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APPARATUS AND METHOD FOR THE IN VIVO DIAGNOSIS OF SKIN 
CONDITIONS BY NEAR INFRARED SPECTROSCOPY 

ABSTRACT 

A novel apparatus and method for collecting and diagnosing skin 
conditions is presented. The apparatus employs a beam of visible and near 
infrared light (400 - 2500 nm), which is passed to the skin via a fiber optic cable. 
Scattered light is collected by the same fiber optic cable, and a spectrum of the 
skin is produced. The method includes collection of spectra from affected skin 
regions (tumors) and non-affected skin regions (normal) in the same patient. 
After taking into account each patient's normal skin spectra, univariate and 
multivariate processing techniques are used to determine differences between 
skin tumors. 

In all cases, normal and tumor spectra were significantly different at 
distinct wavelengths in the spectrum by paired Wests. Classification of near 
infrared spectra of the most difficult clinical diagnoses comparing benign and 
malignant tumors resulted in accuracies between 79.9% and 100%. Clearly this 
NIR spectroscopic technique has the potential for the non-invasive diagnosis of 
skin diseases. 

INTRODUCTION 

Skin cancer is the most common cancer among North American 
populations. In 1999, it is estimated that there will be 70,000 new cases of skin 
cancer in Canada (National Cancer Institute of Canada, 1999) and more than 1 



million new cases in the United States. However, when recognized and removed 
early, skin cancer is curable. Since there often is ambiguity in the clinical 
diagnosis of skin tumors by visual examination, histological analysis of skin 
biopsies is the standard technique for confirmation of a diagnosis. Physicians, 
therefore, must make the difficult decision as to which and how many suspicious 
skin tumors should be biopsied and/or treated. Thus, biopsies * of benign 
neoplasms are often performed in order to err on the side of caution and some 
malignant conditions are overlooked. In addition, long-term monitoring and 
follow-up of tumors with repeated skin biopsies is often required, regardless of 
whether a benign or malignant phenotype is diagnosed. 

A rapid, non-invasive technique that could be utilized for initial screening 
and selection of skin tumors prior to biopsy would be tremendously useful for 
reduction of both patient morbidity and health care costsv Infrared (IR) 
spectroscopy may be that tool. IR spectroscopy uses the translueency of tissues 
and biological materials to IR light, either in the mid* (2500 * 50 000;nm) or near 
(700 - 2500 nm) region, to gain an ^appreciation of tissue biochemistry (Jackson 
et al., 1997). We have shown that biopsies from basal cell carcinoma; squamous 
cell carcinoma and melanocyte tumors have distinct mid IR signatures when 
compared to normal skin (Mcintosh et al, 1999a). Novel processing methods 
have been developed that now make it possible for mid IR spectroscopy to be 
used for the non-subjective diagnosis of skin biopsies (Mcintosh et al., 1999b; 
Mansfield et al., 1999). However, the diagnostic potential of mid IR spectroscopy 
in vivo is limited, as mid IR light has a path-length limitatidrfbf 10-1 5 m . This is 



of little value in vivo, allowing only the very most superficial layers of the stratum 
comeum to be sampled. However, in near IR (NIR) spectroscopy, scattering, 
rather than absorption phenomena, dominate the spectrum. The reduced 
absorption in the NIR therefore allows much greater path-lengths to be utilized, 
offering the possibility of in vivo studies of human tissues. 

Hence, NIR spectroscopy may be the ideal screening methodology for 
skin tumors. With NIR spectroscopy a fiber optic probe is used to illuminate the 
sample; this light penetrates the tissues and scatters throughout, illuminating the 
area under the probe. Light that exits the tissue is collected by the same fiber 
optic probe for analysis. Changes in the amount of light recovered at various 
wavelengths can be indicative of changes in the biochemistry in the illuminated 
region. The NIR measurements are rapid, non-destructive and non-invasive. No 
sample preparation is required, enabling acquisition of spectra directly from the 
skin. 

The NIR region is often sub-divided into the short NIR wavelengths (680- 
1000 nm) and long wavelength NIR (1000-2500 nm) for convenience. At shorter 
NIR wavelengths, the heme proteins (oxyhemoglobin and deoxyhemoglobin) and 
cytochromes dominate the spectra, and their absorptions are indicative of 
regional blood flow and oxygen consumption. Long wavelength NIR absorptions 
arise from overtones and combination bands of the fundamental molecular 
vibrations of C-H, N-H and O-H groups. Therefore, a NIR spectrum gives an 
indication of the distribution and concentration of these molecular species in 
tissue. 



The purpose of this study was to determine whether NIR spectra of 
malignant skin tumors were discemable from benign tumors and normal skin. In 
vivo NIR spectra of six common skin tumors were studied, due to their high 
prevalence and difficult differential diagnosis, namely 1) actinic keratosis, a pre- 
malignant neoplasm of the epidermis. This tumor may develop into the 
malignancy called squamous cell carcinoma. 2) basal cell carcinoma (BCC), a 
malignant skin tumor arising from basal cells of the epidermis. This tumor rarely 
metastasizes, but may cause extensive local growth and spread. 3) dysplastic 
melanocyte tumors and malignant melanoma, which consist of atypical 
J£ melanocytes. Dysplastic melanocyte tumors are usually considered to be pre- 

i - 

\Z malignant neoplasms of the skin, although there is controversy about their 

N 

fn propensity to develop into malignant melanoma. 4) lentigo, a benign skin disease 

Ifl 

\j caused by an increased number of melanocytes. 5) nevi (moles), benign 

5 

M neoplasms of the melanocyte. 6) seborrheic keratosis, a benign,; neoplasm of 

W epidermal cells, and the most common human neoplasm. Univariate and 

111 

K i- multivariate techniques were used to determine* differences between*-the skin 
tumors. 

EXPERIMENTAL 

Subject selection: 

The study base consisted of patients referred for definitive diagnosis of a 
skin tumor to the Dermatology Clinic (JT), and for whom proper patient care 
necessitated a biopsy of theirtumor(s). Upon decistowby the dermatologist that a 
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biopsy(ies) was required, the patient was referred to the study nurse. To be 
included in the study, the subject must have been Caucasian with fair skin and 
18 years of age or older. The subject was excluded from the study if they: 1) 
were using any skin medication on the site of the tumor, 2) were presently 
undergoing radiotherapy or chemotherapy, 3) had either Type I or Type II 
diabetes. Following an explanation and discussion of the Subject Information 
Package (as per NRC's Research Ethics Protocol #1998-14) informed consent 
was obtained. The subject was then enrolled in the study and assigned a subject 
number. 

In total, 153 patients (83 females, 70 males) were enrolled in the study 



\t between the period of October 1998 - June 1999. The mean age of the study 
W population was 59 years (ranging from 21 to 85 years). In total, 195 tumors were 

in 

sampled from the study population of 153 subjects. The mean size of the tumors 



sampled was 8 mm, measured at the tumor's greatest extent. Twenty-five 

O percent (25%) of the neoplasms were located on the trunk, 14% on the upper 

C11 

L -3 limbs, 8% on the lower limbs, and 53% were located on the head and neck. The 

* ** 

•SAX 

subject's skin type was recorded based on the Fitzpatrick Scale of Skin Type 
Classifications (reference). Two percent (2%) of subject's had type I skin (always 
bums easily, never tan), 13% had type II (always bums easily, tan minimally), 
67% had type III (burn moderately, tan normally) and 18% had type IV (bum 
minimally, always tans well). Type V (rarely burns, tan profusely) and Type VI 
(never burns, deeply pigmented.) skin types were not encountered in this study. 
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Acquisition of spectra: 

Near infrared (NIR) spectra were recorded in the 400-2500 nm range in 2 
nm steps using a commercial spectrometer (Foss NIRSystems Model 6500) 
equipped with a bifurcated NIR fiber optic probe from the same*company. The 
system was connected to a laptop computer. A 99% Spectralon reflectance 
standard (LabSphere Inc., North Sutton, NH) was used to convert raw data into 
reflectance spectra. Each reflectance spectrum consisted of 32 co-added scans 
at 10 nm resolution. 

The fiber optic probe was positioned 0.5 mm from the measurement site 
using a Magic Arm (Professional Lighting Supports, Italy) (Figure 1). For all 195 
cases, triplicate NIR spectra were taken from: 1) the tumor and 2) a similar area 
of normal-appearing* skin. Acquisition of each spectrum took 40 ^seconds. Prior to 
obtaining the?readings the subject's skin was&cleansed using a Kimwipe and 70% 
alcohol. The*end*of the probe*was*iCleanedswith»a Kimwipenand^0% alcohol 
between measurement^sites*»and* subjects^Measurements were taken in a 
darkened, quiet room, and the subject was asked4o refrain from talking and to sit 
as still as possible during the measurements. Any movements were recorded. 

After the NIR spectra were acquired, the dermatologist biopsied the 
subject's tumor(s). Biopsies were sent to the pathologist. Slides were evaluated 
by one pathologist (ANN) according to accepted criteria for the histopathological 
diagnosis skin tumors (see Barnhill, 1998 for a complete description of 



histological classification), and detailed reports were issued. The pathological 
diagnosis was used as the gold standard against which to classify the spectra. 



Spectral processing and analysis: 

Prior to any analysis, the spectra were truncated to 400 to 1840 nm, 
leaving a total of 720 data points per spectrum. Triplicate spectra from each site 
were closely examined. Any spectrum that lay outside 2 standard deviations from 
the mean was removed from each group. In all instances "outlier" spectra 
resulted from patient movement, as recorded by the study nurse. The spectra 
W were then averaged by taking the median value for each wavelength. This 
fJJ resulted in one tumor spectrum and one normal spectrum for each case. 

Based on the pathology report, spectra were grouped into one of seven 
categories (see the Introduction for a discussion of the tumor type): 1) actinic 



HI 

m 
in 



m keratosis, 2) basal cell carcinoma, 3) dysplastic melanocyte tumors, 4) lentigo 5) 

a 

O nevi, 6) seborrheic keratosis, 7) other benign tumors. Table 1 summarizes the 
number of cases from each of the above groups. 

Spectra were pre-processed by normalizing each spectral group (tumor 
and normal) using the equation: (original data - sample mean)/standard 
deviation. The normalization routine corrected for offset and minor movement 
variables. Paired t-tests were performed to test any significant differences 
between tumor spectra and their normal counterparts. Subsequently each normal 
spectrum was subtracted from each tumor spectrum in a pair-wise fashion. This 
resulted in one "difference spectrum* for each case, representing tumor only. 
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Linear discriminant analysis (LDA) was performed on difference spectra. 
Prior to LDA, an optimal region selection algorithm, developed by the Informatics 
group at the Institute for Biodiagnostics (Nikulin et al., 1998), was used to select 
the 4-6 most discriminatory regions of the spectra. Optimal regions were used for 
classification by the 'leave-one-out" method of LDA (Eysel et al., 1997, Mansfield 
et al., in press). LDA determines the boundaries that best separate classes by 
computing the distance of each spectrum from all class centroids. It then 
allocates each spectrum to the class whose centroid is nearest. 

RESULTS AND DISCUSSION 

An average NIR spectrum is shown in Figure 2. Spectra are plotted 
showing the amount of light absorbed by the skin at each wavelength between 
400 and 1840 nm. Eachvpeak^in the spectrum can^be^assignedwkto a specific 
substance found in the skin^ The-majon-peak^are assignedtas»shown«in Figure 2. 

A purely qualitative^examination^of ^NIR^speetra«iish0wediKmo«elear^cut 
differences in spectral groups. *Thfe*variability^ from 
different individuals*was*quite^large^due»to*such factors^as^skin^typerage, sun 
exposure, skin location, etc. In particular the inter-person variation was the 
largest in the visible region of the spectrum (400-680 nm), due to distinct 
absorptions from hemoglobin and the pattern of melanin scattering (i.e. skin type 
variation). Therefore, any statistical test performed on skin spectra must take into 
consideration each person's normal NIR spectrum. Since normal spectra were 
taken from a skin-region^ contralateral to that : of the tumor;* we* were able to 
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perform: 1) paired Wests comparing each tumor spectra to its normal counterpart 
(Figure 3) and 2) difference spectra (Figure 4). 

Paired t-tests comparing normal and tumor spectra (by tumor group) 
demonstrated significant differences (p<0.05) at many areas of the spectra, as 
shown in Figure 3. In Figure 3 mean normalized spectra (blue and red traces) are 
shown overlaid on corresponding plots of the p value (black traces) at spectral 
wavelengths. Each tumor type pairing showed a distinct pattern of significance, 
suggesting that success may occur in differentiating the different tumor types. 

However, as demonstrated by the paired t-tests, there are many 
wavelengths of light that are statistically different, and it is difficult to predict 

a 

J* which wavelengths are the most diagnostic. Therefore, the non-subjective 

q multivariate pattern recognition techniques of region selection and LDA were 

Li! 

sj utilized on difference spectra (Figure 4). Difference spectra were obtained by 

M subtracting each tumor-normal pairing for each group shown in Figure 3. Two- 

a 

£3 class LDA's were performed on the clinically difficult classifications (Tables 2-7). 

cn 

^ In LDA tables, the numbers in rows represent the histological classification of the 
tissue, while results in columns represent the calculated classification. The 
bracketed numbers in the LDA tables represent histological classification (rows) 
compared to clinical classifications (columns). 

There are many difficult clinical diagnoses when dealing with skin tumors. 
In this manuscript, only the difficult comparisons between benign and 
premalignant/maiignant tumors are examined. No attempt is made to differentiate 
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between benign lesions or malignant conditions, since in the case of malignant 
conditions, biopsies and removal are ultimately required. 

Differentiating dysplastic melanocyte tumors (malignant) from nevi 
(benign), dysplastic melanocyte tumors (malignant) from lentigo (benign), and 
actinic keratosis (pre-malignant) from lentigo (benign) are some of the- most 
difficult problems in the dermatology clinic. LDA on NIR spectra from these three 
groupings performed excellently, with 91.5% accuracy when differentiating 
dysplastic melanocyte tumors from nevi (Table 2), 100% accuracy when 
differentiating dysplastic melanocyte tumors from lentigo (Table 3) and 84.5% 
accuracy when differentiatng actinic keratosis from lentigo (Table 4). All three of 
these LDAs performed better than the dermatologist's clinincal diagnosis (i.e. 
86.6%, 92.8% and 78.2% respectively). 

Other difficult clinical diagnoses are:* 1) dysplastic melanocyte tumors 

M (malignant) vs. seborrheic keratosis (benign) 2) BCC (malignant) ^seborrheic 

f *^ 

0 keratosis (benign) and 3) actinic keratosis^premalignant) vs. seborrheic^keratosis 
(benign), the results of which are shownnn Tables 5, 6 and 7, respectively. In 
these three comparisons, the dermatologist . performed better-than the LDA (i.e. 
84.7% vs.. 100%, 79.9% vs. 95.8%, 80.3% vs. 95.7%). However, the 
dermatologist still decided to biopsy the tumors, despite his correct clinical 
evaluation. For example, in the instance of dysplastic melanocyte tumors vs. 
seborrheic keratosis (Table 5), the dermatologist biopsied 19 tumors which were 
benign, although his clinical impression was correct. If he had used NIR 



in 
N 

in 
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spectroscopy in combination with LDA, he would have needed to only biopsy 1 
tumor (which is a false positive) out of the original 19. 

Therefore, it is proposed that NIR spectroscopy be used in combination 
with LDA as a method of confirmation of diagnosis, not as a stand-alone 
technique without clinical examination. If used in this manner, the number of 
unnecessary biopsies would be greatly reduced. 

CONCLUSIONS 

We believe this is the first extensive near infrared spectroscopic study of 

- n the most common skin tumors seen in a general dermatology clinic. Clearly this 

S3 

- M NIR spectroscopic technique has the potential for the non-invasive diagnosis of 

in 

£S skin diseases. The results of the present study show that NIR spectroscopy 

\!j differentiates normal skin from that with a tumor. In addition, although 

M preliminary, it appears that NIR spectroscopy may be a promising method for the 

a 

C3 discrimination of malignant from benign skin tumors. 

c 
a 
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Table 1 : Tumor case statistics. 



CASES 


NUMBER 


uauuucu iivjiii aiuuy 




Withdrew 

V V III IUI w V¥ 




No bion^v on NIR ^ite 


12 


iMO iNir^. rriecioijrerrioriio 




Outstanding (awaiting removal) 


26 


Included in study 




Actinic Keratosis 


34 


Ddodl Lrt?ll Odl UK lijl lid 




Dysplastic melanocytic tumors 


10 


Lentigo (actinic and solar) 


13 


Nevi 


22 


Seborrheic Keratosis 


24 


Other benign tumors 


17 




153 


Total 


195 



Table 2: Two class LDA: Dysplastic melanocytic tumors vs. nevi. 





Dysplastic 


Nevi 


% 


specificity 


Positive 


Negative 




melanocyti 




correct 




predictive 


predictiv 




c tumors 








value 


e value 


Dysplastic 


7 


1 


87.5 


95.5 


95.1 


88.4 


melanocyti 


(6) 


(1) 


(85.7) 








c tumors 














Nevi 


1 


21 


95.5 


87.5 


88.4 


95.1 




(2) 


(14) 


(87.5) 









Overall accuracy - 91.5 % (clinical - 86.6 %) 
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Tabl 3: Two class LDA: Dysplastic melanocyte tumors vs. lentigo. 





Dysplastic 


Lentigo 


0/ 
70 


specificity 


rosiiive 


in egauve 




melanocyti 




correct 




predictive 


predictiv 




c tumors 








value 


e value 


Dysplastic 


8 


0 


100 


100 


100 


100 


melanocyti 


(6) 


(1) 


(85.7) 








c tumors 














Lentigo 


0 


12 


100 


100 


100 


100 


(0) 


(6) 


(100) 









Overall accuracy - 100 % (clinical - 92.8 %). 



Table 4: Two class LDA: Actinic keratosis vs. lentigo 





Actinic 


Lentigo 


% 


specificity 


Positive 


Negative 




keratosi 




correct 




predictive 


predictiv 




s 








value 


e value 


Actinic 


31 


2 


93.9 


75.0 


79.0 


92.5 


keratosis 


(28) 


(1) 


(96.5) 








Lentigo 


3 


9 


75.0 


93.9 


92.5 


79.0 


(4) 


(6)' 


(60) 









a Overall accuracy - 84.5 % (clinical «■ 78.2 %) 

M 

Q 



Table 5: Two class LDA: Dysplasticmelanocytie-tumors vs. seborrheic keratosis. 





Dysplastic* 


Seborrhea 


» % 


specificity 


.Positive^* 


Negative 




melanocyte 


c 


correct-- 




predictive 


predictiv 




c tumors 


keratosis 






value 


e value 


Dysplastic 


6 


2 


75.0 


94.4 


93.1 


79.1 


melanocyti 


(6) 


(0) 


(100) 








c tumors 














Seborrheic 


1 


17 


94.4 


75.0 


79.1 


93.1 


Keratosis 


(0) 


(19) 


(100) 









Overall accuracy - 84.7 % (clinical - 100 %) 
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Table 6: Two class LDA: Basal cell carcinoma (BCC) vs. seborrheic keratosis. 





BCC 


Seborrhei 
c 

keratosis 


% 

correct 


specificity 


Positive 

predictive 

value 


Negative 
predictiv 
e value 


BCC 


27 
(29) 


2 

(1) 


93.1 
(96.7) 


66.7 


73.6 


90.6 


Seborrheic 
keratosis 


6 

d) 


12 
(19) 


66.7 
(95) 


93.1 


90.6 


73.6 


Overall accuracy - 79.9 % (clinical - 


95.8 %) 








Table 7: Two class LDA: Actinic keratosis vs. seborrheic keratosis. 






Actinic 

keratosi 

s 


Seborrhei 
c 

keratosis 


% 

correct 


specificity 


Positive 

predictive 

value 


Negative 
predictiv 
e value 


Actinic 
keratosis 


31 
(28) 


2 

(1) 


93.9 
(96.5) 


66.7 


73.8 


91.7 


Seborrheic 
keratosis 


6 

(1) 


12 

(19) 


66.7 
(95) 


93.9 


91.7 


73.8 



Overall accuracy - 80.3 % (clinical - 95.7 %) 



DESCRIPTION OF FIGURES 



Figure 1 : A schematic of the experimental set-up. A light source emits near 
infrared light, which is directed to the skin's surface via fiber optics. Reflected 
light is collected via fiber optics and directed to the detector, which subsequently 
is processed to produce a spectrum. 



Figure 2: A mean (+/- SD) normal near infrared skin spectrum. The major peaks 
are assigned. 
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Figur 3: Paired West results comparing normal and skin tumor spectra. The 
mean normalized spectra (blue and red traces) are shown overlaid on 
corresponding plots of the p value (black traces) at spectral wavelengths 
between 400 and 1840 nm. 

Figure 4: Difference spectra from skin tumors. Difference spectra were obtained 
by subtracting each tumor-normal pairing for each group shown in Figure 3. 
Actinic keratosis (blue), basal cell carcinoma (red), lentigo (green), dysplastic 
nevi (black), nevi (pink) and seborrheic keratosis (brown) are shown. 

We Claim: 

1 . A method for the diagnosis of skin conditions, comprising exposing an 
area of a patient's skin to IR light, collecting the scattered light, producing an IR 
spectrum of the scattered light, comparing the*spe,etrum of scattered light to a 
spectrum of IR light collected from another skin»area from the same patient, and 
determining the differences between the skin areas* 

2. An apparatus for the in vivo diagnosis of skin conditions, comprising light 
source means for producing a beam of IR light, means for directing the light to 
the skin surface of a patient and for collecting the scattered light to produce an IR 
spectrum thereof, and computer means for analyzing the spectrum. 

3. Apparatus for taking measurements of diseases of the skin using near 
infrared light, comprising: 

a light source for producing a beam of infrared light, 
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a fiber optic for directing the light to the skin's surface and for collecting 
the reflected light, thus producing an infrared spectrum, and 
a computer for analyzing the produced spectrum. 



4. A method of diagnosing diseases of the skin using near infrared light with 

the apparatus in claim 1, which includes: 

collection of baseline normal spectra from each patient from an area of 

non-affected normal skin, 

collection of abnormal spectra from each patient's skin disease, 
providing a database of spectra previously collected from normal and 

abnormal skin, including placing spectra in the database in groups based upon a 

thorough histological examination of skin biopsies which were taken after in vivo 
JJj near infrared spectroscopy was performed, and 

}.i comparison of the patient spectra to the database spectra and 

v assigning the patient spectra to a group, thereby diagnosing the skin condition 

m 
in 

sg 5. The method according to claim 4, wherein the patient spectra are 
f . collected by placing a fiber optic probe above the skin. 

fl 

rn 6. The method according to claim 4, wherein comparison of the patient 
^ spectra to the database spectra takes into account each patient's normal skin 
type. 

7. The method according to claim 4 wherein the patient spectrum is 
compared to the database spectra over wavelengths of high accuracy, the 
wavelengths of high accuracy being selected from the group consisting of 
wavelengths corresponding to oxyhemoglobin, deoxyhemoglobin, water, proteins 
and lipids 

8. The method according to claim 4 wherein the patient spectrum is a mean 
spectrum of two or more spectra from each normal and each abnormal skin site. 
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